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We report on a systematic treatment of Coulomb corrections for 3- and n-particle 
(Bose-Einstein) correlations, leading to a removal of a 100 % systematic error from 
earlier Coulomb corrections of 5-particle Bose-Einstein correlations in heavy ion 
collisions. These results are generalized to include strong final state interactions 
and core- halo effects. 

1 Introduction 

One of the most important tasks of high energy heavy-ion studies is to prove 
the existence of the elusive quark-gluon plasma and to study the proper- 
ties of this predicted new state of matterl. Hanbury-Brown Twiss (HBT) 
interferometryB of identical particles has become an important tool as it 
can Jbe rUspd to measure the evolving geometry of the interaction region, see 
refsBu'ErQi] for some of the recent reviews of this rapidly expanding field. The 
quantitative interpretation of the HBT-results depends critically on the under- 
standing of the role of the Coulomb interaction between the selected particles. 
However, Coulomb interactions between three (or more) charged particles are 
notoriously difficult to handle, as evidenced by the many decades of research 
in atomic physics aimed at obtaining solutions of the three-body Coulomb 
scattering problem which are accurate and possess a wide range of applicabil- 
ity. 

Starting from a three-body Coulomb wave function which is at least ex- 
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act for sufficiently large inter-particle separations, we determined the effect of 
Coulomb final state interactions on the three-particle Bose-Einstein correla- 
tion function of similarly charged particles in ref.u. We estimated numerically 
that the familiar Riverside approximation is not precise enough to determine 
the three-body Coulomb correction factor in the correlation function if the 
characteristic HBT radius parameter is between 5 and 10 fm, which is the 
range of interest in high-energy heavy ion physics. 

We then generalized the method oLCoulomb wave function corrections 
from three to n charged particles in ref.Eil Here, we re-formulate this work so 
as to be able to include the possible effects of i) additional strong final state 
interactions and ii) core-halo model corrections. 



2 Coulomb and final state interaction effects in ?i-particle 
correlations 

The n-particle Bose-Einstein correlation function is defined as 

r (v , v _ iv„(ki,---,k») 

where N n (ki, • • • , k„) is the n-particle, and iVi(k^) the single-particle inclusive 
invariant momentum distribution. The three-momentum vector of particle i 
is denoted by k^. ft is quite remarkable that this complicated object, that 
carries quantum mechanical information on the phase-space distribution of 
particle production as well as on possible partial coherence of the source, 
can be expressed in a relatively simple, straight-forward manner both in the 
analytically solvable pion-laser model of refs.EJOEjLas well as in the generic 
boosted-current formalism of Gyulassy and Padulao as 

n 

]rn G ( k " k -.) 
c " (ki ' ' ' ' ' kn) = nluG^) ' (2) 

where cr( n ' stands for the set of permutations of indices (1, 2, ■ ■ • , n) and <Ji 
denotes that element which replaces element i in a given permutation from 
the set of u'"' . Regardless of the details of the two different derivations 

G(k J ,k,) = (a t (k 4 )a(k,)) (3) 

stands for the expectation value of a^(ki)a(kj). 
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In the relativistic Wigner- function formalism, in the plane wave approxi- 
mation G(ki,k2) can be rewritten as 



G(ki,k 2 ) = J d*xS(x,Ki 2 ) exp(% 2 ■ x) (4) 
i^i2 = 0.5(fci + fc 2 ), qi2 = k 1 -k 2 , (5) 



where a four- vector notation is introduced, fej = ( \J + k? , kj ) , and a ■ b 
stands for the inner product of four-vectors. Due to the mass-shell constraints 
— W to? + hf , G depends only on 6 independent momentum components. 
In any given frame, the boost-invariant decomposition of Eq. (Q) can be rewrit- 
ten in the following, seemingly not invariant form: 

G(ki,k 2 ) = J d 3 x S Kl2 (x) exp(-iqi 2 x), (6) 

Sk u W = / ' dtexp(ip Kl3 q 12 t) S(x,t,K 12 ), (7) 
/3 ifl2 = (k 1 +k 2 )/( J B 1+J B 2 ). (8) 

If n particles are emitted with similar momenta so that their n-particle Bose- 
Einstein correlation functions may be non-trivial, Eqs. (^,||) will form the basis 
for evaluation of the Coulomb and strong final state interaction effects on the 
observables. On this level, all the correlations are build up from correlations 
of pairs of particles. 

In order to treat the Coulomb corrections to the n-particle correlation 
function exactly, knowledge of the n-body Coulomb scattering wave function is 
required. We restrict ourselves to the case that the transverse momenta of all 
n particles in the final state in their common center of mass are small enough 
to make a nonrelativistic approach sensible. Hence the problem consists in 
finding the solution of the n-charged particle Schrodinger equation when all 
n particles are in the continuum. 

Consider n distinguishable particles with masses rrii and charges e,, % = 
1,2, • • •, n. Let x$ denote the coordinate (three-)vector of particle i. Then 
Yij = Xj — Xj is the relative coordinate between particles i and j, and k^ = 
(mjkj — TOjkj) / (rrii + rrij) the canonically conjugate relative momentum. 

The n-particle Schrodinger equation reads as 



#o+ E ^-S^Lt ) ..k„(xi,---,x„) = 0, (9) 

«i=i 
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where 

" k 2 

E = T.^>° < 10 » 

is the total kinetic energy for n particles in the continuum. Hq is the free 
Hamilton operator and 

V i j(rij) = V?(r ij ) + V?(r ij ) (11) 

the interaction between particles i and j, consisting of a short-range strong 
potential V$ and the long-range Coulomb interaction, V^(r.y) = e.%e.jl \ Yij |. 

Although for n = 3 an exact numerical solution of the Schrikhnger equa- 
tion @ with = V^ 7 has been achieved recently for E > 0E_3, because of 
its complexitiy the resulting wave function can not easily be used for gen- 
eral purposes. In addition, for n > 3 a numerical solutiop-.is beyond present 
means. For a brief discussion of the related difficulties seetZl. This means that 
for practical purposes analytic, even if only approximate, n-charged particle 
wave functions are still needed. 

For n = 3 such wave functions are available in the form of the explicit solu- 
tions of the Schrodiagec. equation in all asymptotic regions of the three-particle 
configuration spacalllij. The simplest form, which is also easily generalized 
to arbitrary particle numbers, applies to the (dominant) asymptotic region 
conventionally denoted by f2o and characterized by the fact that p. roughly 
speaking - all three inter-particle distances become uniformly largetj. In the 
final states of heavy-ion reactions, where a large number of charged particle 
tracks appear, the mutual, macroscopically large separation of tracks is one 
of the criteria of a clean measurement. This suggests that in order to study 
Coulomb effects on n-body correlation functions, knowledge of the wave func- 
tion in f^Q 11 ^ , the region in n-particle configuration space where all interparticle 
distances become uniformly large, i.e., | — > oo for all values of (ij), may be 
sufficient. Hence, an appropriate, approximate n-particle Coulomb scattering 
wave function is 

n 

<>.. )kB (*i, •••,*») - n (i2) 

t<j=i 

Here, tp k \ r ij) is the continuum solution of the two-body Coulomb 
Schrodinger equation (with a reduced mass of fiij = mirrij j {rn% + rrij)): 

2 
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Figure 1. Coulomb wave function correction (filled circles) as compared to the less sub- 
stantiated Gamow correction (squares) for 5-particle Bose-Einstein correlations for a source 
sizes R = 1 fm (left panel) and R = 5 fm (right panel). The inset shows the ratio of these 
two correction factors. 



describing the relative motion of the two particles i and j with energy kf -/2/Zjj. 
The explicit solution is 

= % <' k /' "/,,• l;i(| k*j II T tj | -kyr^)], (14) 

with Nij = e~ 7rr,ii / 2 T(l + and T]ij — etej^ij/ | | being the Coulomb 

parameter. F[a, b; x] is the confluent hypergeometric function and r(x) the 
Gamma function. The foundation of Eq. ( [l2| ) in fundamental scattering theory 
and its expected range of validity are discussed in ref.EI 

Additional strong final state interactions can be taken into account in 
a straight-forward manner by substituting in Eq. ( |l2| ) for the two-particle 
Coulomb wave functions tp k . (r^) the solutions of the Schrodinger equation 
with the Coulomb plus strong potential, 

In this way a few dominant, strong interaction induced phase-shifts can be 
taken into account for the pair under consideration. On the level of the two- 
parlicle correlation function, a similar technique has already been applied in 
reffi 
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3 Application to high-energy heavy-ion and particle collisions 



The correlation function measuring the enhanced probability for emission of 
n identical Bose particles is given by Eq. (1). This correlation function is 
usually, due to meager statistics, only measured as a function of the Lorentz 
invariant Q n , defined by the relation 

a 

Ql = E ( fc *-%) 2 - ( 16 ) 

i<j=l 

We can now calculate the Coulomb effects on the n-particle correlation 
function using 

xnr=i^(xi)|< )s k „( X i,---,x, 

Kcoulomb(Qn) — J - (^) 



Here, p(x ; ) is the density distribution of the source for particle i (normalized 
to the total number of particles), taken as a Gaussian distribution of width R 

in all three spatial directions and VP^..^ (xi, • • • , x„) ~ n"<j=i e tkijrij is the 
n-body wave function without any final state interaction. The superscript S 
indicates appropriate symmetrisation for identical particles. This formulation 
makes it possible to extract information on the source size R, and to compare 
this value with that extracted by means of a generalized n-particle Gamow 
approximation via K ( ^ ulomb {Q n ) = Yli<j^i^ij> where G ^o =1 N ij I 2 - 

Pion n-tuples were sampled randomljoE2l from tha.NA44 data sample of 
three pion events produced in S-Pb collisions at CERNlEl. Fig. 1 demonstrates 
that in case of a characteristic effective source size of 1 fm, the difference 
between the n-particle Gamow and the Coulomb wave function corrections 
are smaller than 10 % for n = 5 particles. However, with increasing num- 
ber of particles and/or with increasing effective source sizes, this difference 
increases dramatically. For instance, for 5 particles the naive generalized 
Gamow method overestimates the Coulomb correction as compared to the 
much better substantiated Coulomb wave function integration method by a 
factor of 2. 



4 Core-halo corrections 

The results given in the previous section have to be corrected for the ef- 
fects that may arise from the existence of a halo of long-lived hadronic 
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resonancesc3S. In this case, the effective source function p(x) has two com- 
ponents, one pertinent to a smaller region corresponding to the core of the 
interactions and one to a larger region that corresponds to the (unresolvable) 
halo: 

p(x) = p c (x) + p h (x) (18) 
(n) = (n c ) + (n h ) (19) 

(rii) = dxpi(x), i = c,h. (20) 



The effective intercept parameter of the two-particle Bose-Einstein corre- 
lation function reads as 

A* = {(n c )/{n)f (21) 

if the core fraction is independent of the relative momentum. If the core-halo 
model is applicable to a given data set, the radius Rh of the halo does not 
matter as long as it is big enough. For instance, it can be a Gaussian with a 
radius of 20-40 fm; in any case Rh >> h/Qmin, the inverse of the two-particle 
relative momentum resolution. 



5 Summary 

In this contribution we reviewed the state-of-art method for treating Coulomb 
corrections to 3- and n-body correlation functions. For small effective source 
sizes of about 1 fm, the pairwise product of Gamow factors gives a good ap- 
proximation for the n-body Coulomb correction factor. However, for source 
sizes typical in high energy heavy ion experiments, a more substantiated cal- 
culation based on the pairwise product of relative Coulomb wave functions 
is required in order to remove a possible 100 % systematic error from this 
correction. 

We have also generalized this correction method for the inclusion of strong 
final state interactions and pointed out that the Coulomb corrections have to 
be performed self-consistently with the core-halo corrections; in particular 
the effective intercept parameter of the two-particle Bose-Einstein correla- 
tion function has to be taken into account when evaluating the 3- or n-body 
Coulomb corrections. 
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